Association of a healthy beverage score with total mortality in the adult population of Spain: A nationwide cohort study

Background Despite the substantial evidence of the relationship between diet and mortality, the role of beverage consumption patterns is not well known. The aim of this study was to assess the association of the adherence to a Healthy Beverage Score (HBS) and all-cause mortality in a representative sample of the Spanish adult population. Methods and findings We conducted an observational cohort study using data from the Study on Nutrition and Cardiovascular Risk in Spain (ENRICA), which included 12,161 community-dwelling individuals aged ≥18 years recruited in 2008 to 2010 and followed until January 2022. At baseline, food consumption was collected using a validated diet history. The HBS consists of 7 items, each of which is scored from 1 to 4 (highest adherence). The HBS ranges from 7 to 28 points with a higher score representing a healthier pattern. Adherence was assigned as a higher consumption of low-fat milk, and coffee and tea, a lower consumption of whole-fat milk, no consumption of fruit juice, artificially sweetened beverages, or sugar-sweetened beverages, and no or moderate consumption of alcohol. Total mortality was ascertained by linkage to the Spanish National Death Index. Statistical analyses were performed with Cox models and adjusted for the main confounders, including sociodemographic, lifestyle, dietary variables, and morbidity. After a mean follow-up of 12.5 years (SD: 1.7; range: 0.5 to 12.9), a total of 967 deaths occurred. For all-cause mortality, the fully adjusted hazard ratio (HR) for the highest versus lowest sex-specific quartiles of HBS was 0.72 (95% confidence interval [0.57, 0.91], p linear-trend = 0.015), corresponding to an 8.3% reduction in the absolute risk of death. A linear relationship between the risk of death and the adherence to the HBS was observed using restricted cubic splines. The results were robust to sensitivity analyses. The main limitation was that repeated measurements on beverage consumption were not available and beverage consumption could have changed during follow-up. Conclusions In this study, we observed that higher adherence to the HBS was associated with lower total mortality. Adherence to a healthy beverage pattern could play a role in the prevention of premature mortality.


What did the researchers do and find?
• We included a representative sample of 12,161 adults (18 years and older) from Spain who were recruited in 2008 to 2010 and followed up until 2022.A total of 967 deaths occurred.
• Participants were categorized according to their adherence to the HBS.
• A higher total score was achieved with a higher consumption of low-fat milk, and coffee and tea, no consumption of whole-fat milk, fruit juice, artificially sweetened beverages, sugar-sweetened beverages, and no consumption or moderate consumption of alcohol.
• Each HBS item scored from 1 (minimum adherence) to 4 points (maximum adherence) and the HBS ranged from 7 to 28 points.The higher the HBS, the healthier.
• When comparing extreme categories, higher adherence to the HBS was associated with lower all-cause mortality in the Spanish adult population, with an 8.3% reduction in the absolute risk of death.

What do these findings mean?
• The adherence to the HBS could serve as a potential diet-based strategy to prevent premature mortality.

Introduction
The influence of unhealthy dietary factors on adverse outcomes, including premature mortality, is a public health concern [1].Thus, the association between diet and mortality has been examined by using several approaches such as analyzing individual nutrients, food and, more recently, assessing dietary patterns and indexes [2][3][4].Most of the indexes include mainly solid food (e.g., meat, poultry, fish, as well as fruit and vegetables) [5], although some of them also comprise beverages [6].
The mechanisms by which beverages influence health are complex and are not only based on the nutritional quality of their components (such as energy provided, macronutrients, fiber, minerals, and vitamins) [7], but also rely on other factors such as satiety mechanisms and factors that affect the assimilation of beverages such as rapid gastric emptying and intestinal absorption [8].Moreover, the addition of artificial sweeteners influences mortality [9].On the other hand, some beverages are also an important source of other additives (e.g., phosphates) as well as contaminants from packaging or processing (e.g., organophosphate esters, phthalates) [10].
In 2015, a healthy beverage index, based on commonly consumed drinks, was developed to evaluate the role of beverage quality on cardiometabolic risk in adult Americans.A low adherence to this index was associated with several detrimental cardiometabolic markers [11].Consistent results were obtained in another US study where the association between the adherence to a healthy beverage pattern and total mortality was evaluated.Data were obtained from a cohort of 2,283 adults, aged �21 years, with a previous diagnosis of mild to moderate chronic renal insufficiency.A healthier beverage index was inversely associated with the progression of chronic kidney disease and all-cause mortality [12].
Although the association of specific beverages has been studied previously, the role of a healthy beverage index and its association with mortality has not been assessed in the general population yet.We hypothesized that higher adherence to a 7-item Healthy Beverage Score (HBS), previously proposed by Hu and colleagues [12] and adapted to the Spanish beverage consumption, could be associated with lower mortality.Therefore, the aim of this study is to assess the association between the HBS and all-cause mortality in a representative cohort of Spanish adults.

Study design and participants
Data were obtained from the Study on Nutrition and Cardiovascular Risk in Spain (ENRICA) whose methods have been reported elsewhere [13].In brief, 13,105 individuals aged �18 years were recruited from 2008 to 2010.A stratified cluster sampling based on the census sections of Spain was performed to guarantee the representativeness of the sample.Sample weights were based on the size of municipalities, sex, and age.Three sequential stages were followed for data collection.First, sociodemographic, lifestyle characteristics, and morbidity information was obtained through a telephone interview.Second, blood and urine samples were collected on a first home visit.Third, a physical examination and a face-to-face dietary history (DH-EN-RICA) were completed during a second home visit.The response rate was 51% and the main reasons for non-participation were refusal to provide a blood sample (51.7%) and not being interested in the study (37.8%).
The Clinical Research Ethics Committee of La Paz University Hospital in Madrid provided ethical approval.All participants from the ENRICA Study gave written informed consent after explaining the details of the study.

Study variables
Dietary history.Information on diet was obtained through a computerized dietary history (DH-ENRICA), conducted by trained and certified nonmedical interviewers.The DH-ENRICA collected information on 861 items of food, with 82 beverages included.Participants informed about all items of food and beverages consumed at least once every 2 weeks in the previous year.Food consumed during weekdays and weekends were considered.A total of 127 sets of digitalized photos, household measurements, as well as the usual proportion sizes of food from typical Spanish recipes were used to estimate portion sizes in grams per day.Regarding beverages, a total of 14 digitized photos and 23 household measurements were used to later estimate beverage consumption in milliliters per day.In addition, for alcoholic beverages, the consumption of ethanol in grams per day was calculated using Spanish food composition tables [14].The validity correlation coefficients in HD-EN-RICA for beverages were: 0.71 for coffee, 0.69 for milk, 0.40 for soft drinks, and 0.64 for alcoholic beverages [15].
The Healthy Beverage Score (HBS).A Healthy Beverage Score (HBS) was previously described by Hu and colleagues [12].Based on the HBS, we built a 7-item HBS modifying its cut-off points to fit with the beverage consumption of a representative sample of the Spanish adult population.Each item of the HBS scored from 1 (minimal adherence) to 4 points (maximal adherence) based on sex-specific categories of consumption.Thus, the HBS ranged from 7 (low adherence) to 28 points (high adherence).The higher the HBS, the healthier the pattern.Items were grouped in 2 main components: adequacy and moderation (Table 1).Two beverages were considered as adequacy components: low-fat milk as well as coffee and tea consumption.For these 2 components, the higher the score, the healthier the pattern.No low-fat milk consumption scored 1, while the remaining sample was divided into tertiles among consumers; coffee and tea consumption was grouped into quartiles.Five items were included as moderation components: whole-fat milk, fruit juice, artificially sweetened beverages, sugarsweetened beverages, and alcohol.For these 5 items, the higher the consumption, the lower the score, with a specific classification for alcohol consumption.Whole-fat milk and sugar-sweetened beverages scored 4 for no consumption and the remaining sample was divided into tertiles among consumers.Fruit juice and artificially sweetened beverages consumption scored 4 for no consumption and 1 for any consumption.The scoring of fruit juices and artificially sweetened beverages was decided on the basis of the lack of a wide range of consumption, and to maintain the relative weight of these items in the score, as previously described by Hu and colleagues [12].Finally, for alcohol consumption, participants with no consumption or moderate drinkers (<40 g/day for males and <24 g/day for females) scored 4, and heavy drinkers (�40 g/day for males and �24 g/day for females) scored 1.
Mortality assessment.All-cause mortality was ascertained through a computerized linkage with the Spanish National Death Index.Participants were followed from baseline in 2008 to 2010 to January 31, 2022.Follow-up was censored at the date of death or at the end of follow-up, whichever occurred first.
Confounders.Participants provided information regarding age, sex, educational level, and smoking status which were obtained through the computer-assisted telephone interview performed at baseline.On the second home visit and following standardized procedures, blood pressure, and weight and height were measured.Body mass index (BMI) was calculated as weight divided by the square of the height in meters (kg/m 2 ).Leisure time and household physical activity were evaluated with the EPIC short questionnaire, collecting information on 17 activities.Then, each activity was multiplied by their respective energy expenditure rate in metabolic equivalents in hours per week (METs-hour/week) [16] and total energy expenditure was obtained by summing up all activities.Hours spent watching television were used to account for sedentary activities.In order to control for good dietary quality, analyses were adjusted for total energy intake, fiber, fruit and vegetable consumption, as well as the Mediterranean Diet Adherence defined by Trichopoulou and colleagues [17] without including alcohol.Blood samples collected on the first home visit were centrally analyzed in the CORE laboratory of La Paz University Hospital in Madrid.A colorimetric enzymatic method with lipase and glycerol kinase (for triglycerides) and a colorimetric enzymatic method with cholesterol-oxidase, esterase, and peroxidase (for cholesterol) were used.To define hypertriglyceridemia, we used a threshold of �150 mg/d in fasting plasma triglycerides levels, and for hypercholesterolemia, a fasting plasma total cholesterol level of �200 mg/dL or prescribed lipid-lowering medications.Hypertension was defined as �140/90 mmHg or taking antihypertensive medication.Analyses were also controlled for the number of chronic conditions (chronic obstructive pulmonary disease, coronary heart disease, stroke, heart failure, osteoarthritis, cancer, depression diagnosed by a physician, and diabetes), as well as the number of prescribed medications to consider prevalent morbidity.

Statistical analysis
Across sex-specific quartiles of adherence to the HBS, age-adjusted baseline characteristics of participants were computed using marginals.Age-adjusted means for continuous variables and age-adjusted proportions (%) for categorical variables were provided.To estimate hazard ratios (HR) and their 95% confidence intervals (CIs), Cox proportional hazards regression models were built and age was considered as the underlying time metric.The survey command was applied to account for the complex sampling design.The lowest category of adherence to the HBS was used as reference.Three sequential Cox models were used.Model 1 was an unadjusted model.Model 2 was additionally adjusted for age, sex, educational level, smoking status, ex-drinker status, BMI, physical activity in leisure time, total energy intake, fruit and vegetable consumption, total fiber intake, hypertriglyceridemia, hypercholesterolemia, hypertension, number of selfreported chronic conditions, and number of medications.Finally, Model 3 was adjusted for the Mediterranean index by Trichopoulou excluding alcohol (maximum score = 8), but excluding fruit, vegetable, and fiber consumption.Schoenfeld residuals were plotted against time to assess proportional hazards assumptions and visually no violations were found.To test for linear trend, categories of the HBS were modeled as a continuous variable.The doseresponse relationship was assessed by using restricted cubic splines with 3 knots (at the 10th, 50th, and 90th percentiles).Interactions between the HBS and age (<65 years versus �65 years), sex (male versus female), BMI (<30 kg/m2 versus �30 kg/m2), physical activity (�median 61.5 METs-hour/week versus >median 61.5 METs-hour/week), vegetable consumption (�median 183.5 g/d versus >median 183.5 g/d), adherence to the Mediterranean diet without including alcohol (�median 4 versus >median 4), and the prevalence of chronic conditions (yes/no) were tested by including multiplicative terms in Model 3. Sensitivity analyses were conducted excluding deaths in the first 3 years of follow-up to account for the effect of subclinical conditions at baseline.Also, individual items of the HBS were assessed according to Model 3 plus adjustment for the remaining items that are part of the score.
This was a preplanned study and data analysis was conducted according to a prespecified plan (S1 Text).The HBS was previously used to assess the association between adherence to HBS and age-related frailty in a sample of Spanish older adults [19].A minor modification was introduced in alcohol consumption classification.For older adults, moderate alcohol consumption was considered the healthy option due to their high cardiovascular risk [19].However, in the current analysis, which involves the adult general population (aged �18 years) with lower cardiovascular risk, the category of alcohol considered healthy was "no consumption or moderate consumption." This study was reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline (S1 STROBE Checklist).
A two-tailed p value less than 0.05 was considered statistically significant.All analyses were performed with Stata, version 17.0 (StataCorp, College Station, Texas, United States of America).

Results
The median age of participants (N = 12,161) was 46 years old (interquartile range 35 to 61) and 52.6% were females.Compared with those in quartile 1 (less healthy) of adherence to the HBS, participants in quartile 4 (healthier) were older, more frequently females, with a higher level of education and with a less sedentary lifestyle, and were more physically active.Also, those in quartile 4 had lower energy intake, consumed more fiber, fruit and vegetables, showed a higher adherence to the Mediterranean diet, and had more frequently hypercholesterolemia (Table 2).Sex-specific cut-off points for individual items of the HBS are shown in S1 Table .In accordance with the rules for the construction of the HBS, compared with those in quartile 1 (less healthy), participants in quartile 4 (healthier) consumed more low-fat milk, coffee and tea, and alcohol, but consumed less whole-fat milk, fruit juice, artificially sweetened beverages, and sugar-sweetened beverages (Table 3).After a mean follow-up of 12.5 years (SD: 1.7; range: 0.5 to 12.9) and 151,459 person-years of follow-up, a total of 967 deaths occurred.The HR for all-cause mortality when comparing extreme quartiles of the adherence to the HBS was 0.72 (95% CI, 0.57 to 0.91, p for linear trend = 0.015) in the fully adjusted model (Table 4).The decrease in absolute risk of death was 4.3% for quartile 2, 6.3% for quartile 3, and 8.3% for quartile 4. No significant interactions were found for age, sex, BMI, physical activity, vegetable consumption, or adherence to the Mediterranean diet (without including alcohol).However, a statistically significant interaction was found when stratifying for the presence of at least 1 chronic condition (p = 0.030).Among those with at least 1 chronic condition, higher adherence to the HBS was associated with lower mortality.No association was observed among those with no chronic conditions (Table 5).After excluding the first 3 years of follow-up, the inverse association between the adherence to the HBS and total mortality remained similar (S2 Table ).When assessing dose-response, a linear relationship was observed using restricted cubic splines (p value for non-linearity = 0.010) (Fig 1).When individual HBS items were analyzed using Model 3 for adjustment, a higher consumption of coffee and tea, and no consumption of fruit juices and artificially sweetened beverages contributed most to the association with lower all-cause mortality (

Discussion
In this large population-based study of Spanish adults with a mean follow-up of 12.5 years, a higher adherence to the HBS was inversely associated with total mortality, after adjusting for potential confounders.Those with higher adherence to the HBS had an 8.3% reduction in the absolute risk of death compared to those with lower adherence.The association was linear and robust.It may also be of particular interest to people with preexisting chronic conditions, as they had lower mortality, although these findings need to be confirmed in future research.
Regarding to items of the HBS, 2 recent prospective studies performed in Spain found an inverse association between coffee consumption and all-cause mortality [20,21].Results from the EPIC study (with 500,000 participants from 10 European countries) [22] and from the UK Biobank study were also similar to the findings in this study [23].Our results are also in line with meta-analyses comprising cross-sectional studies, longitudinal cohorts, as well as interventional studies [24][25][26].The beneficial effect of coffee might rely, among others, on the antioxidant and anti-inflammatory activity exhibited by its bioactive components, mainly melanoidins, chlorogenic acids, and caffeine [27].These compounds reduce oxidative stress and inflammation [28], enhance endothelial function [29], and counteract carcinogenesis on in vitro studies [30].Coffee also increases the metabolic rate [31], improves the glucose metabolism [32], and lowers long-term blood pressure [33].Moreover, coffee could reduce mortality even in those with impaired caffeine metabolism [34] and independently to the addition of sweeteners [35].However, high coffee consumption has been associated with an increase in serum levels of total cholesterol, LDL-cholesterol, and triglycerides [36].On the other hand, evidence suggests that coffee consumption above 4 cups/day is not associated with further lower mortality [25].Spain is included among the European countries with the lowest tea consumption [37] and we did not find studies that evaluated its relationship with mortality among Spanish adults.Therefore, it is unlikely that tea consumption accounts for our results.However, in literature, both all-cause and cardiovascular mortality were reduced among tea consumers [38].
The effect of milk on health has been widely studied due to its fatty acid composition [39].Two recent cohort studies showed that low-fat milk consumption was associated with lower all-cause mortality when compared to whole-fat milk consumption [40,41].Whole-fat milk has a higher content of saturated fats that has been related to an increase in LDL-cholesterol and atherosclerosis [42].As a result, whole-fat milk consumption could be particularly harmful among individuals with known cardiovascular risk.However, in a clinical trial among normocholesterolemic individuals, whole-fat milk consumption showed no impairment in lipid profile nor in glucose-insulin metabolism when compared to low-fat milk consumption [43].It is of note that, milk could modulate satiety mechanisms [44] and also has several components with potential beneficial effects, such as caseins with antioxidant properties [45].Milk is also rich in minerals, vitamins, and other bioactive compounds involved in anti-inflammatory and immune regulation [46].A recent meta-analysis showed that whole-fat milk consumption was associated with a higher risk of all-cause, cardiovascular disease and cancer mortality; however, low-fat milk showed a protective but nonsignificant association [47].
Results on fruit juice consumption and mortality mostly depend on the distinction between processed or fresh fruit juice [48].A recent meta-analysis showed that processed fruit juice consumption was associated with a higher risk of type 2 diabetes and total mortality [49].Evidence on fresh fruit juice consumption and health, however, is insufficient to draw conclusions.Results of a cohort study from the US with 13,440 participants showed that a higher consumption of 100% fruit juice was associated with a higher mortality [50].Conversely, a study with 198,285 individuals from the UK found a positive association between sugar-sweetened beverages and mortality, but not for 100% fruit juice consumption [51].Similarly, a recent meta-analysis of prospective cohorts concluded that there was no association between 100% fruit juice consumption and all-cause mortality [52].Several studies have also found ).Adjusted as in Model 3. Data were adjusted for age (years, continuous), sex (male, female), educational level (primary or less, secondary, university), smoking (non-smoker, former smoker, current smoker), ex-drinker (yes/no), BMI (<25, �25 and �30, >30 kg/ m 2 ), time watching TV (hours, continuous), physical activity (METs-hour/week, continuous), energy intake (kcal/day, continuous), hypertriglyceridemia (yes/ no), hypercholesterolemia (yes/no), hypertension (yes/no), number of chronic conditions (0, 1, and �2), number of medications (0, 1-3, >3), adherence to the Mediterranean diet without including alcohol (maximum score = 8), and for the rest of items of the HBS (as appropriate).Age was the underlying time metric.CI, confidence interval; HBS, Healthy Beverage Score; HR, hazard ratio.https://doi.org/10.1371/journal.pmed.1004337.g002that, compared with 100% fruit juice, sugary or processed fruit beverages produce harmful glucose levels after ingestion, mainly due to the higher content of free sugars [53].In our study, bottled, sweetened, as well as fresh fruit juices were analyzed together as a unique item because of their rapid absorption [8] and similar effect on postprandial glucose levels [54].Fructose intake, particularly from sugar-sweetened beverages at any dose, or from fruit juice at higher doses, contributes a rapid extra dietary energy source that could explain its detrimental effect on health [55].However, food-based dietary guidelines from various countries from Europe Union, including Spain, consent to replace occasionally 1 daily portion of fruit with fresh fruit juice [56,57].
In order to lower calorie intake and control body weight, artificially sweetened beverages could be adequate short-term substitutes [58].However, when considering the long-term influence of artificially sweetened beverage consumption, several studies have found associations with higher obesity, hypertension, type 2 diabetes, stroke, cardiovascular disease incidence and mortality, and all-cause mortality [59][60][61].Since these beverages contain few to no calories nor sugars [51], some investigations have related them with weight gain as a result of an increased consumption of sweet food due to a greater affinity for sweet flavors or the perception of eating fewer calories [62].In addition, their flavoring components have been associated with the formation of advanced glycation end-products [63], which are involved in the development of metabolic diseases [64].Moreover, some sweeteners such as sucralose and saccharin could induce glucose intolerance and alterations in gut microbiota [65] that are linked to obesity [66].
In literature, a low to moderate alcohol consumption is related to a reduction in all-cause mortality [67].Biological explanations for this protective role on health are based on lipid regulation, insulin response, and endothelial function [68] resulting from the modulation of some anti-inflammatory biomarkers [69].However, at high doses, alcohol is detrimental to cardiovascular health and is related to several types of cancer [70].A harmful alcohol consumption is associated to neurodegenerative processes [71], microbial dysbiosis [70], and an increased intestinal permeability that leads to a permanent hepatic exposure to bacterial translocation, oxidative stress, and other inflammatory components [72].In addition, alcohol use could result in hepatic steatosis and de novo lipogenesis, and also could reduce the utilization of lipids [73].Lastly, alcohol may injure myocardium with potential cardiomyopathy and heart failure [74], and increase the risk of hypertension [75].On the other hand, the beneficial association of alcohol consumption with mortality found in some studies may rely on abstinence bias, insufficient adjustment for covariates, or consumption changes due to disease detection [76,77].However, recent studies from an epigenetic perspective have proposed that alcohol at restricted doses could be particularly beneficial in older adults based on changes in alcohol metabolism related to age [78].In this regard, a meta-analysis for the Global Burden of Disease Study has proposed a change from sex-specific to age-specific recommendations on alcohol consumption [79].Then, low alcohol consumption could be beneficial among older adults, but not for younger adults.For our analyses, we considered that being a heavy drinker was harmful because of its well-established association.
In a previous study, a 10-item Healthy Beverage Index was constructed using an a priori approach based on the US recommendations for beverage consumption [11].Similar to this index, Hu and colleagues described the HBS as a more suitable score for use in large epidemiological studies [12].The HBS excluded water consumption as well as 2 items on total energy from beverages and calculations of daily fluid intake.In our study, the same scoring weights (from 1 to 4) were maintained for all items in the HBS.However, in contrast to Hu and colleagues, we considered both no alcohol consumption and moderate alcohol consumption as healthy.Additionally, we also modified the HBS cut-off points of the items to fit with the beverage consumption of the Spanish adult population.As a result, the HBS used in our study retained the same items as originally described by Hu and colleagues, as well as the relative weight of the items.
The use of the HBS as an overall measure of beverage consumption has several advantages.First, the use of this score overcomes the limitations of analyses of relationships between individual beverages and diseases, as beverage consumption may be correlated, and an increase in consumption of one beverage may be associated with a decrease in the others.Secondly, the HBS could be a complementary tool for assessing adherence to dietary patterns that include only solid foods, in order to assess dietary quality as a whole.Thus, the HBS could serve as a simple and rapid screener to obtain information on the quality of beverage consumption from the general population, similar to other indexes used to assess adherence to certain diets, such as the Mediterranean diet.Also, the use of this 7-item pattern may be an optimal choice when dealing with patients in time-constrained clinical settings.Finally, the HBS includes items on commonly consumed beverages and could be easily adapted to other populations with only minor modifications to account for their specific beverage consumption.
We have used the HBS in the general population and caution should be exercised in deriving beverage consumption recommendations from this score in specific populations, especially those with restricted fluid requirements, long-term liquid diets, and other preexisting conditions involving fluid consumption.Further studies are also needed in specific population subgroups.In addition, a future study could consider intercorrelations and specific populationbased patterns of beverage consumption using an a posteriori approach.
This study has some limitations.First, when measuring diet, non-differential misclassification of the exposure is always possible, in general, resulting in an underestimation of the associations found.Second, no information concerning behavioral changes or repeated measurements on beverage consumption were available, and beverage consumption could have changed during follow-up.Third, water consumption was not available in this study.Water is universally recommended as a safe beverage and as the main source of hydration.As water does not provide energy, macronutrients or micronutrients, its consumption is considered free for the general population.
There are also some strengths.To our knowledge, this is the first examination on the relation between a healthy beverage score and all-cause mortality among the Spanish adult population.In addition, we used a dietary history that allowed us to collect information on beverages with validity and reproducibility in a Spanish population.Also, the national vital statistics records, accessed through linkage to the Spanish National Death Index, ensured an extensive follow-up of the cohort for mortality assessment.Finally, several confounders were considered in more adjusted models.
In conclusion, in this representative study of the Spanish adult population, higher adherence to the HBS was associated with a reduction in total mortality.As the consumption of a healthy solid diet should be encouraged, adherence to a healthy beverage consumption pattern may also play an important role in the prevention of premature mortality as part of public health nutrition prevention strategies.

Table 1 .
Scoring criteria for the HBS in the ENRICA Study (2008 Fig 2).Unadjusted results are also shown (S3 Fig).

Table 2 . Age-adjusted baseline characteristics of participants in the ENRICA Study (2008-2010) by quartiles of the HBS (N = 12,161). Characteristics Quartiles of adherence to the HBS p for trend
Chronic conditions included: chronic obstructive pulmonary disease, coronary heart disease, stroke, heart failure, osteoarthritis, cancer, depression diagnosed by a physician, and diabetes.HBS, Healthy Beverage Score; BMI, body mass index; METs-hour/week, metabolic equivalents in hours per week. a

Table 3 . Beverage consumption by quartiles of the HBS in the ENRICA Study (2008-2010) (N = 12,161). HBS components Quartiles of adherence to the HBS a p value
bAlcohol was defined as the consumption of ethanol in grams per day.HBS, Healthy Beverage Score; SD, standard deviation.https://doi.org/10.1371/journal.pmed.1004337.t003